Two s-wave solutions for superconductivity in the extended Hubbard model 
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The existence of more than one solution for s-wave pairing in the extended Hubbard model is not 
often realized. This possibility was noted by Friedberg et al. [Phys. Rev. B50, 10190 (1994)] in 
the case of two electrons on a lattice, without further analysis. In the present paper, the second 
solution is found also in the case of superconductivity of extended s-wave symmetry in the extended 
Hubbard model. The properties of both s-wave solutions are examined by mean-field methods and 
thresholds for their appearance are given. A possibility of the first-order transition is discovered and 
modifications of the phase diagram are calculated. Results of this paper in the limit of low electron 
density are also applicable to the bound, two-electron pairs. 
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I. INTRODUCTION 



The discovery of high-temperature superconductivity 
(HTS) has triggered an intensive effort to reveal its 
mechanisms. This aim has not been attained yet but 
the amount of knowledge concerning HTS is substantial. 
From this wealth of data certain key issues emerge, which 
must be addressed by the eventual microscopic theory. 
One of them is question of the symmetry of the order 
parameter. Instead of originally accepted s-wave, nowa- 
days it is believed that the hole-doped HTS materials 
are predominantly of the dx'^-y^-wave symmetryjiiSi Nev- 
ertheless the research concerning the s-wave symmetry is 
not to be abandoned. The s-wave symmetry is found in 
fuUerideSf^ perovskites^ and MgB2^ not from the HTS 
group but very interesting materials though. The or- 
der parameter of s symmetry can appear as subdomi- 
nant one in d-wave superconductors due to orthorombic 
distortions of the crystal&L^ There is evidence of the s 
type HTS in electron-doped materialsSiiSiiiiiS and pos- 
sibility of the s to d change of the symmetry of the or- 
der parameter with increasing doping in the hole-doped 
materials >i^*i^ Another scenario with s-wave symmetry 
includes the possibility of different order parameters in 
the bulk and on the surface of HTS superconductorii^ 
As there are reports of different mechanisms responsi- 
ble for pairing and phase coherence there is possibility of 
pseudogap having s symmetry, different from the d sym- 
metry of the order parameter liSiii Of other possibilities 
we can not exclude possibility of finding an s-wave sym- 
metry pairing in heavy fermions, which are notorious for 
multiplicity of phases. 

There are few rigorous results concerning s-wave su- 
perconductivity or BCS in general. We know that in 
the dilute limit BCS mean-field equations go over to 
the Schrodinger equation for the bound two-electron 
pair,-^'^^ i.e. in this limit the BCS equations become 
exact Another result was given by Randeria et al}^ 
stating that creating two-electron bound state of s-wave 
symmetry is a necessary and sufficient condition for BCS- 
type superconductivity to exist in two dimensions in a 
dilute limit. The results of applying mean-field approxi- 



mations to the Hubbard-type Hamiltonians as well as the 
properties of the bound pairs are well known (e.g. Refs 

13,20^21,22^23,24^25-^^ Q^^^, q£ ^J^g gg^Qm UOticcd aSpCCtS of 

mean-field treatment of the extended s-wave pairing and 
bound pairs of this symmetry is subject of this paper. 



II. FORMALISM 

Let's consider extended Hubbard model in the stan- 
dard notation: 

ija i ijaa' 

(1) 

where we sum over nearest neighbors (n.n.), U and W are 
on- and intersite interactions respectively, tij is nearest- 
neighbor hopping integral and /x - chemical potential. In 
the mean-field approach the above Hamiltonian takes the 



fori 



JX. 



H = 2J(£'£ ^ t^)nka - 2j(^fec[-fcLj,j^ + H.c.) + const. 



ka 



(2) 

where, using order parameter < c-qiCq^ >, we intro- 
duced a gap Afc! 



^ = 77) ]Vkq<C-qlCq^ > 
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(3) 


Vkq = -U - W-tk-q 


(4) 


7fc = ) ]exp{ik-S) 


(5) 


£k = -tlk 


(6) 


JI = fi- {U/2 + Wz)n 


(7) 



z is coordination number; for the full discussion of the 
model the reader is referred to the RefiS& It is easy to 
derive the self-consistent equation for the energy gap: 



^^ = ^E^^-9^tanh/3£;,/2 



(8) 



where the quasiparticle energy: 



E„ 



y^(e,-7i)2 + |A,| 



W 



(9) 



and (3 — 1/ [ksT) where T is temperature and ks Boltz- 
mann constant. To solve the model we usually make the 
ansatztii 



Afe = Ao + A^7a; + A^T^fc 



(10) 



where rjk = 2(cos/ea; — cosfcj^) and particular terms re- 
fer to on-site s-, extended s- and d-wave pairings. Self- 
consistent equations for the gap separate into the s-wave 
and d-wave part. 

Finally what is left for solving for general, non-zero U 
and W in the ground state is a set of three self-consistent 
equations for Aq, A-y and chemical potential /I: 
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A. = --^E7.(^o + 7,A.)^ 
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in the case of the extended s-wave pairing, and 
~1^ ^^''2E'n 



(11) 
(12) 

(13) 
(14) 



in the case of the d-wave pairing. In the following we 
will focus on the s-wave pairing only. The results de- 
scribed in the remainder of the paper are calculated for 
the rectangular density of states (DOS) with z = 4, so 
they approximate two dimensional square lattice. For 
the reasons that will be explained later on the author 
supposes that they are of more general nature. 



III. RESULTS 



Equations (|ll() - (|13|l in the dilute limit yield 
Schrodinger equation for a bound pair, as was told 
in the Introduction. As is well known, in two dimensions 
infinitesimally small on-site attraction creates bound 
pair for ly = 0. For W^ > a threshold must be crossed 
to obtain a bound state. Analytic formula for the 
threshold and U-W phase diagram are given for example 
jjjji3j27 rpj^g j^ggg known facts are as follows: 

(i) The threshold formula for the gap, analogical to the 
one given in,^'^-^^ can be derived from BCS-MFA equa- 
tions (|11|I - H13() without resorting to the limit n ^ 0. Us- 
ing rectangular density of states it reads: 



Wr. 
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At 



l + {n-lY{l + lQt/U) 



(15) 
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FIG. 1: Critical values for existence of extended s-wave su- 
perconductivity calculated for the rectangular DOS, for elec- 
tron densities n = (also n = 2, full line) and n = 0.2 (also 
n = 1.8, dashed line). Dot-dashed and full line with circles 
are thresholds for n = 0.2, calculated for the exact DOS in 
two dimensions (2d) and in three dimensions (3d, simple cu- 
bic lattice) respectively. Exact 2d threshold for n = is the 
same as for rectangular DOS. Exact 3d threshold for n = 
is given by solid line with squares. Dotted lines are the axes 
of coordination system. Dot-dashed lines with triangles are 
crossover boundaries for n = 0.2 and exact 2d DOS. Crossover 
boundaries for n = coincide with the threshold lines in 2d. 
Full, flat, horizontal line is a threshold for the d-wave pair- 
ing, calculated for n = 0. Numbers "0", "1" and "2" are the 
names of the areas described in the text. U and W are in the 
half-bandwidth units D = 4f in 2d and D = 6f in 3d. 



Let's note that the limit n ^ yields correct result of 
RefsiiiSi The use of the exact density of states does not 
bring changes to the value of Wcrit for ti — > but the 
changes appear and increase with growing n, as shown in 
Fig. [2 

(ii) The formula (|15|l is valid on the entire U-W plane, 
both for positive and negative values of U and W . The 
cited papers show the figure of the fourth quarter of the 
coordination system only. The full plot is given in Fig. ^ 
As we can see in the third quarter of the coordination 
system, for U and W attractive enough, a second thresh- 
old line is obtained! (In the case of two-electron bound 
pair second solution was also noted in Ref.'^^ but without 
further analysis). 

This way the area in Fig. ^ can be divided into three 
main parts: Area "0" - to the right of the rightmost 
threshold line, mostly first quarter (C/ > 0, M^ > 0) and 
small parts of the second and fourth quarter of the coordi- 
nation system (not exceeding asymptotic values: Uas/t — 
-16(n-l)7(l + (n-l)2) and W^s/t - -4/(l + (n-l)2)). 
There are no bound states and no gap for these param- 
eters, the solutions describe antibonding states. Area 
" 1" ~ between both threshold lines, mostly second and 
fourth quarter - in this area there is only one solution of 



the equations (|ll|l - l|13|l . which exists even for infinite on- 
or intersite repulsion. The further from the rightmost 
threshold line we are the stronger is the binding between 
electrons in the pairs. Points of the same binding en- 
ergy form a line approximately parallel to the rightmost 
threshold line. Area "2" - below the second threshold 
line, in the third quarter - for parameters from this area 
two solutions of equations (|ll|) - (|13|l exist: one of the 
type of the "area 1", with large binding energy, given 
by a "distance" from the first threshold line; the other, 
with smaller binding energy and smaller gap, given by 
the "distance" from the second threshold line. We will 
call these solutions type 1 and type 2, respectively. Let's 
note, that the range of U, for which both types of pairing 
exist, increase with increasing n (the asymptote moves 
from Uas/4:t = —2 for n = to Uas = for n = 1) while 
analogical range of W decrease (the asymptote moves 
from Was/'^t = -1/2 for n = to Was /it = -1 for 
n = 1). The properties of both types of pairing are de- 
scribed in the following. 

For fixed W, for U -^ — oo, type 1 solution displays 
what we will call here "pure-s" type behavior: |Ao| ~ U, 
A^^O,-p^ U/2 (if H/ > then A^ < 0, if W" < then 
A-y > 0). When we move on the U axis from left to right, 
towards smaller values of jt/j then, for certain intermedi- 
ate value of |[/| (not too large if W is not large and at- 
tractive) this "pure-s" behavior disappears. If 14^ > Was 
then with increasing U we must reach threshold value 
Ucrit and for U -^ Ucrit Aq decreases non-linear ly to 0, 
I A-yl goes through extremum and also decreases to while 
fl goes to constant, non-zero value. liW < Was and so- 
lution of type 1 exists in the whole range of U from — cx) 



to -hoo, then for U 



-cxD all three quantities Aq, A- 



and JL go to the constant values: Aq to a small negative 
one and A-y to a larger, positive one (see Fig. O. When 
M^ < and \W\ is large enough, then "pure-s" behavior 
disappears ior U ~ W instead of small to intermediate 
values of |C/|, as shown in Fig. [3 

The behavior of type 2 solution is somewhat comple- 
mentary: for U — > — cx), when type 1 solution shows 
"pure-s" behavior, type 2 solution is constant. Moving 
right on the U axis type 2 solution goes through a transi- 
tion for the same range of U as type 1 solution, and from 
now on it behaves in "pure-s" way until U reaches the 
threshold for type 2 solutions, where Aq and A^, disap- 
pear and JI = const. Let's note, that the slope of linear 
dependence of Aq vs U differs just by the sign from the 
analogical slope of type 1 solution and that A.y is of op- 
posite sign to Aq in type 2 solution. 

The constant values reached by the type 1 solution for 
[/ — > -l-oo are the same as the ones reached by the type 2 
solution for U -^ — oo. These constants for small n and 
\W\ large enough can be calculated as: 

jL = W/2 (16) 

Ay/Ao--|7l| (17) 

Aq « -0.6V^ (18) 

As shown in Fig.|2|for fixed U, when W is varied, instead 
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FIG. 2: Type 1 (full line, full line with circles) and type 
2 (dashed hne, dashed line with circles) solutions of the ex- 
tended Hubbard model vs [//4f for W/U = -30 and n = 0.01. 
Upper panel shows Aq and A^ and lower panel Jl. Dotted 
lines show axes of the coordination system. In the inset the 
enlargement of the area near the origin of the coordination 
system in the upper panel is shown. 



of "pure-s" behavior type 1 solution displays "pure ex- 



tended s" behavior for W 



-oo: 



\Wl Ao ^ 0, 



-p r^ W/2. Apart from that, the behavior of solutions vs 
W is analogous to the behavior of the solutions vs U: for 
W -^ -l-oo type 1 solution goes to the constant values 
(now A-y to the small negative one while Aq to the big- 
ger, positive one) and for W -^ — oo type 2 solution goes 
to the same constants. For small n and \U\ large enough 
these constants are: 



M 



Ao/A, 



U/2 



(19) 
(20) 
(21) 



Large values of fixed parameters in Figs 13 and |21 were 
chosen to emphasize the characteristic behavior of the 
examined quantities. The plot of chemical potential Jl vs 
U for W/At — —30 is practically indistinguishable from 
the plot of 7^ vs W^ for U/At = —30 in the scale of the 
picture (apart from the different critical values for which 
Jl in the type 2 solutions appears) so it was not shown 
again in Fig. 13 

A point to notice is that the gap parameter Aq remains 
finite even in the limit f/ — » oo, despite the fact that on- 
site pairing amplitude < cLcL > must vanish in this 
UmitaSiSa 

Let's note that the considerations above answer the un- 
spoken question, why the equations for extended s-wave 
pairing yield formally (and numerically) the same solu- 
tion in the two very different limits: infinite attraction 
and infinite on- or intersite repulsion. Despite they are 




FIG. 3: Type 1 (full line, full line with circles) and type 
2 (dashed line, dashed line with circles) solutions of the ex- 
tended Hubbard model vs W/U for U/U = -30 and n = 0.01. 
Dotted lines show axes of the coordination system. The de- 
pendence of /I on W is indistinguishable in the scale of the 
picture from the dependence oijl on U in Fig. |21 (apart from 
different thresholds for type 2 solutions), so it was not shown. 
In the inset the enlargement of the area near the origin of the 
coordination system is shown. 



numerically the same, nevertheless they are two different 
solutions of two different types of pairing! 

In the limit of large |[/| or large \W\ the solutions for 
Aq and A.^ are symmetric or antisymmetric regarding to 
the change n ^ 2 — n. For large, negative U (or W) and 
fixed W (or U) Aq (A-y) of "pure-s" (of "pure extended 
s") type 1 solution is symmetric, reaching maximum for 
the half- filled band, while A^ (Aq) is antisymmetric: 
positive (or negative) for less than half-filled band and 
of opposite sign for more than half-filled band, crossing 
through for n = 1. In the opposite limit U (or W) going 
to -l-oo the behavior of Aq and A^ in "type 1" solutions 
change to the reverse: Aq (A^) becomes antisymmetric 
and A^ (Aq) symmetric. The behavior vs n of type 2 
solutions for [/ or H^ ^ — oo is the same as behavior of 
respective "pure-s" or "pure extended s" type 1 solutions 
in the limit U or W ^ +oo. This is the correct behavior 
if the limiting, constant values reached by the type 1 so- 
lutions are to be identical to the type 2 limiting, constant 
solutions for every n. For smaller values of the parame- 
ters \U\, \W\ there is an area in parameter space of "type 
2" solution, where both Aq and A^ are symmetric in the 
band, though of opposite sign. 

A peculiarity of solutions is the existence of critical 
charge density Uc- For a given value of the fixed parame- 
ter {U or W) there exists Uc such that for n in the range 
{ric, 2 — Uc) the transition between the "pure (extended) 
s" part of one solution and the "constant" part of the 
same solution is no longer smooth, but both parts join in 
an non-continuous way: we have a first-order transition. 



where A's and derivative of JI jump. Some approximate 
values of nc{U/4:t,W/4:t) for type 1 solution, given as a 
function of parameters for which the jump takes place, 
are: nc(-29.1, -30) « 0.0475, nc(-9.1,-10) « 0.15, 
nc(-2.75,-3.75) « 0.5, nc{-l.l,-2) « 0.8. We can see 
here that for \W\ large enough and attractive the jump 
takes place for U « W, as was stated before. The values 
of |t/| and \W\ required to obtain any given ric are much 
larger in type 2 solution than in type 1 solution. 

To complete the examination of the n-dependence let's 
notice that at half-filled band solutions fully separate 
and we obtain true pure-s solution, with A^ — and 
no threshold for pairing and true pure extended s solu- 
tion with Ao = and threshold W = -At {W « -3.24i 
for exact 2d density of states) iS& 

Another point to notice is that we can add the third 
and fourth solution to the two found, by changing Aq 
and A-y to — Ap and — A^. The number of nontrivially- 
different solutions seems to depend on the number of cou- 
pled terms from the ansatz Eq. (fTn|l and does not depend 
on the lattice dimensionality. In Fig. 1 there is a plot of 
threshold curves for simple cubic, three dimensional (3d) 
lattice. Their behavior is analogous to the 2d case. In 
particular the threshold line for the type-1 pairing goes 
through the point {U, W) = (0, 0), even for small electron 
densities (except from n = 0). This does not contradict 
the fact that ior U = we have critical value for |W^| to 
bound electron pair in 3d. The existence of bound pairs 
is prerequisite for superconductivity in a dilute limit only 
in 2d^-. A curve for two-electron bound pair in 3d can be 
obtained from Eqs. Hll(l - p3|l . by solving them for n — 0. 
This curve is also given in Fig. 1, plotted in a certain 
distance from other 3d curves for n > 0, in agreement 
with Refsii21 

There appears a question: which of the two solutions 
will be physically realized? The type 1 pairing has larger 
binding energy, higher critical temperature and lower free 
energy so this is the solution to be found even in the area 
of parameter space where both solutions exist. In this 
area care should be taken not to mix the two solutions 
especially for U « W, where the absolute values of solu- 
tions of type 1 become close to the respective values of 
solutions of type 2 and also in the limit of large negative 
U or W, where equations can easily yield type 2 solution 
numerically. 

Wc must also remember that only the relative stabil- 
ity of two extended s-wave solutions was examined in the 
above considerations. The question of stability of another 
symmetry pairings, e.g. d- or p-wave, was not addressed 
at all. While p-wave pairing can be diminished by adding 
antiferromagnetic coupling to the Hamiltonian, d-wave 
can not. From Fig. 1 we can see that for large enough 
\W\, d-wave pairing can surely have bigger binding en- 
ergy than second extended s-wave solution. As we know 
it can win the competition even with the first extended 
s-wave solution (with the lower free energy - see e.gsM). 
In establishing the detailed phase diagram we must thus 
take into account the competition of pairings with other 



symmetries, including s -\- d or s + id. 

Another question is how physical is the first order tran- 
sition found, isn't it just an artefact of the MFA, espe- 
cially that it appears for larger n? This question could 
be probably answered by numerical simulations. If this 
effect is physical it should have influence on the all su- 
perconducting properties of the superconductor with ex- 
tended s-wave superconductivity. In particular one can 
ask whether the found first order transition could be an 
argument in favor of non-smooth crossover between BCS 
and local pairs limits in the models of extended s-wave 
superconductivity. It seems unlikely as we notice that 
\U\ and \W\ being the arguments of ric cited before in the 
text, decrease for growing Uc while the BCS-Bose bound- 
ary (calculated form the Leggett's criterion?^ i.e., that 
chemical potential crosses the bottom of the band) moves 
toward larger values of |C/| and \W\ with increasing elec- 
tron density, so that the two boundaries do not coincide. 
Crossover boundaries for n — and n — 0.2 are plotted 
in Fig. n] There are two branches of crossover bound- 
ary corresponding to the two threshold lines. For n = 
they coincide with the threshold lines, with increasing n 
they move left, towards larger and larger values of |C/| 
and \W\. Again, if other than s-wave symmetry pairing 
is realized, this problem is irrelevant. 

The third issue concerns the possible phase separation 
for large, values of \U\ or \W\. This question was ad- 
dressed by Hartree type analysis e.g. in the paper of 
Robaszkiewicz and Pawlowski^'^, from where the phase 
boundaries shown in Fig. 0] come (solid lines). Under the 
lowest, horizontal line, for large n.n. attraction, we have 
a phase separation of normal state and electron droplets. 
In the area between the solid lines there is s-wave super- 
conductivity (on-site only, denoted as SS) and above the 
upper solid line there is inhomogeneous phase separation 
of charge density waves (CDW) and superconductivity 
(CDW/SS). For the half filled band the ground state is 
CDW ior W > while CDW and superconductivity are 
degenerated for W — 0. 

Our second superconducting s-wave solution does not 
enter this picture but for intermediate values of on-site 
attraction, Eq. ((T^ helps us to make an upper bound- 
ary for superconducting state, denoted by dashed lines 
in Fig. 21 Above this line no superconductivity exists. 
This way for small electron densities, for W large enough 
(above the critical value) a normal state is obtained (de- 
noted as NO in the figure). With increasing n we can 
have a transition to SS and then to CDW/SS state. 
For certain larger W we can have NO-CDW/SS transi- 
tion and for still larger W we obtain NO-CDW-CDW/SS 
transition. This way a phase diagram of Ref.^^ is sub- 
stantially modified. Considering full extended s-wave 
superconducting solution, with order parameter depend- 
ing on both on- and intersite interaction, changes phase 
boundaries in Fig. 0] only in negligible way. 



Let's note, that the described modification of the phase 
diagram takes place only for intermediate values of \U\. 
For larger values of \U\ or more precisely for U < Uas, 




FIG. 4: The phase diagram of the extended Hubbard model 
solved in the mean-field approximation for f//4f — —0.8 
(lines with no symbols) and C//4t = —0.4 (lines with cir- 
cles). Full lines are taken from the paper of Robaszkiewicz 
and Pawlowskii-^ Dashed lines come from Eq. 1151 . N/ED 
denote phase separation between normal state and electron 
droplets, SS - superconductivity, CDW - charge density waves, 
CDW/SS - inhomogeneous phase separation between CDW 
and SS, NO - normal state. 



where Uas was defined in connection with Fig. 1, s-wave 
superconductivity exists independently of W, and there 
is no normal phase area in the phase diagram. 

In conclusion a mean-field criterion for existence of s- 
wave superconducting solution in the extended Hubbard 
model was calculated. The modified W — n phase di- 
agram was shown. A second s-wave solution was found 
and a threshold for its existence together with a full phase 
diagram in the UW plane were calculated. The second 
solution helped to answer the question why limits U (or 
W) = ±oo are formally the same. (These results apply 
also to the two-electron bound states of the same sym- 
metry). First-order transition in the extended s-wave 
solution was found and its possible implications on the 
crossover problem were discussed. The results may apply 
to all properties of the superconducting state of extended 
s-wave symmetry. 
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